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■ Abstract. We discuss the phase structure and fluctuations of conserved charges in 

I two flavor QCD. The importance of the density fluctuations to probe the existence of 

I ' the critical end point is summarized. The role of these fluctuations to identify the first 

. order phase transition in the presence of spinodal phase separation is also discussed. 
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^ ■ 1. Introduction 
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CN I One of the essential predictions of QCD is the existence of a phase boundary in the 



(T, /iq)-plane that separates the chirally broken hadronic phase from chirally symmetric 



O 

O I quark-gluon plasma phase [Ij. Arguments based on effective model calculations [1-12] 
I indicate that at large /i^ the transition along this boundary line can be first order. In 
I I the opposite limit of high temperature and low baryon number density the transition 
from hadronic to quark gluon plasma phase is expected to be continuous or second order 
^ ' depending on the number of quark flavors and the quark masses [7-12]. This suggests 
that the QCD phase diagram can exhibit a critical endpoint where the line of first order 
^ I phase transitions matches that of second order or analytical crossover [13]. The critical 
I properties of this second order critical endpoint are expected to be determined by the 
Ising model universality class [5[ [71 [12]. The existence of a critical end point in QCD 
has recently been studied in lattice calculations at non-vanishing chemical potential in 
2- and (2+l)-fiavor QCD [H [ig [I6] . 

In this paper we discuss how a generic QCD phase diagram appears in effective 
chiral model calculations. We consider the properties of the fluctuations of the conserved 
charges to identify the position of the phase boundary and the critical end point. We 
compare the model results with recent LGT findings in 2-flavor QCD and discuss their 
physical interpretation. Finally, considering deviations from an idealized equilibrium 
picture of the first order chiral phase transition we discuss the properties of charge 
density fluctuations by including spinodal instabilities. 

2. Charge density fluctuations in effective chiral models and in LGT 



To study thermodynamics and the phase structure related with the chiral symmetry 
in QCD we adopt the Nambu-Jona-Lasinio model formulated in the mean field 
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Figure 1. The chiral condensate M (left-hand figure) and the chiral susceptibihty 
Xch (right-hand figure) calculated in the NJL model in the chiral limit as a function 
of temperature for different values of fiq — 300, 275 and 200 MeV corresponding to a 
1st order, the tricritical point and a 2nd order transition [S] . 



approximation ^ i9j . In the chiral hmit this model describes the chiral phase transition 
from a phase of massive quarks with dynamically generated mass M(T, Hq) dependent 
on temperature T and quark chemical potential /Xg to a phase of massless quarks. 

The effective quark mass M(T, ^q) is a measure of the thermal expectation value of 
the quark condensate at finite T and /ig. Thus, M(T, /i^) acts as an order parameter for 
chiral symmetry restoration. In the chirally broken phase M 7^ and it vanishes when 
the symmetry is restored. In the Fig. 1-left we show the T-dependence of the dynamical 
quark mass for different values of /Xg. For small iiq the condensate continuously melts 
in the narrow temperature range indicating the second order phase transition. At large 
the condensate M drops discontinuously from a finite to zero value indicating the 
first order nature of the chiral transition. The fluctuations of the order parameter 
Xch =< > — < M >^, shown in the Fig. 1-right, are finite when crossing the 
first order transition whereas they diverge at the second order line. The divergence is 
directly linked to the appearance of a massless mode in the scalar-isoscalar, so called 
sigma channel, if the transition is second order. The resulting phase diagram obtained 
in the NJL model with the particular set of the model parameters and in the chiral 
limit is shown in the Fig.2-left. It exhibits the generic, QCD-like structure discussed in 
the introduction. For a finite quark mass the second order transition is converted to a 
cross-over transition [1]. 

The existence of the tricritical point (TCP) is seen in Fig. 2-left to be strictly 
related with the strength of the four-fermion interactions in the isovector-vector channel 
controlled by the coupling Gy in the NJL Lagrangian. For large Gy > O.QGs (with Gs 
describing the strength of the scalar four-fermion interactions) the TCP disappears from 
the phase diagram due to strong repulsive interactions in the medium. 

The position of the TCP can be identified by considering the fluctuations of the 
charge density Uq that are quantified by the corresponding susceptibilitie Xq = duq/dfiq 
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Figure 2. The left-hand figure: The NJL model phase diagram in the chiral limit for 
different values of the isovector- vector coupling Gy. The right-hand figure shows the 
net quark number fluctuations for jiq — 310,305 and 300 MeV [S]- 




Figure 3. The left-hand figure: the isovector fluctuations versus T for fixed /iq chosen 
as in the Fig. 2. The right-hand figure: the quark number susceptibility Xq near TCP 
versus t = \T — Tc\/Tc when approaching the TCP along the line of the 1st and the 
2nd order transition as well as at constant T . Also shown is Xq when crossing the 0(4) 
critical line at fixed fiq [9]. 



[9l[T8l[T9]. Fig. 2-right shows the T-dependence of the net quark number susceptibihties 
in the vicinity and at the TCP. The Xq varies rapidly with T and drops discontinuously 
at the phase boundary. The size of this discontinuity grows with increasing /ig up to 
the TCP where the susceptibihty diverges. Beyond the TCP the discontinuity is again 
finite. At /ig = the discontinuity vanishes and the susceptibihty shows a weaker, 
non-analytic structure at the transition temperature. Also at finite quark mass the 
discontinuity along the cross-over line and at the critical end point is replaced by a peak 
which diverges at the critical end point. The critical properties of Xq the TCP and 
along the 2nd order line are universal and are governed by different critical exponents. 
Fig. 3-right illustrates the critical behavior of Xq near the 2nd order critical line 
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Figure 4. The left-hand figure: the isovector xi and the quark number Xq 
susceptibilities at = 0. The right-hand figure: the fourth-order cumulant moments. 
The resuhs are from LGT calculations with A^/ = 2 and ~ 80 MeV [HIIT]. The 
lines for T/To < 1 are the hadron resonance gas model results. 



and at the TCP. In the mean field approximation the critical dependence of Xq on the 
reduced temperature t = \T — Tc\/Tc is consistent with that obtained in the Landau 
theory [71 [9]. For paths approaching the TCP asymptotically tangential to the phase 
boundary, the susceptibility diverges with the correlation length critical exponent 7^ = 1. 
For other paths the critical exponent is 7^ = 1/2. Along the second order critical line, 
the susceptibility remains finite. Going beyond the mean field and including quantum 
fiuctuations in the lagrangian the above critical exponents are renormalized to those 
belonging to the universality class of the 3D Ising model [18]. On the other hand, the 
finite value of Xq along the 2nd order line is consistent with the 0(4) universality class 
[3 [T71 [IHl [19]. There, the critical properties of Xq are governed by the specific heat 
critical exponents a ~ —0.2 fT, T7]||| Consequently, at /ig = the T-dependence of 
Xq is only due to the regular part of the free energy [TJ [T7|. At finite fig the singular 
contribution to Xq results in the cusp structure the strength of which increases with 
fj,q. The fourth-order cumulant d^ ~ d'^P/dfig develops a cusp already at /i^ = and 
diverges at finite fig along the 0(4)-line [T7] . 

Contrary to Xg, the isovector fiuctuations xi shown in Fig. 3-left are neither singular 
nor discontinuous at the chiral phase transition. The non-singular behavior of xi at 
the TCP is consistent with the observation that there is no mixing between isovector 
excitations and the isoscalar sigma field due to the SU(2)y isospin symmetry [19]. 

Recently, the net quark as well as isovector fiuctuations and higher moments were 
computed [El [17] in two fiavor lattice QCD (see Figs. 4 and 5). At fig = 0, Xg shows a 
strong suppression of the fiuctuations in confined phase and a cusp-like structure that 
increases with fig in the vicinity of the transition temperature at fig 7^ 0. Also observed 
on the lattice was a cusp in the forth-order cumulant moment at = 0. The lattice 

I In the NJL model we obtain the mean-field value for this critical exponent, a — 0. 
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Figure 5. The net quark number Xq (left-hand figure) and isovector xi (right- 
hand figure) susceptibilities for different values of p^qjT . The results are from LGT 
calculations with Nf = 2 and rriq ~ 80 MeV [M] . 



results confirmed that the quark fluctuations in the isovector channel xii contrary to 
Xg, do not show any peak structure and depend only weakly on /ig. These properties 
of Xq ^'iid xi ^^^G consistent with 0(4) universality arguments and should be there when 
approaching the critical end point with increasing ^q. However, the above behavior of 
Xq and xi can be also interpreted in terms of the regular part of the free energy due to 
the enhanced contribution of resonances in the vicinity of the transition temperature. 
The LGT results on T and /i^ dependence of Xq and xi as well as c/^ can in the confined 
phase be successfully described [IT] by the resonance gas partition function as seen in 
Figs. 4. This is also illustrated in Fig. 6, where the ratio 2 = d\/d'^q of the fourth- to 
the second-order cumulants as well as the ratio of the quark density to the susceptibility 
are seen to be consistent with the hadron resonance gas model. The observed change 
in cumulants at is due to deconfinement resulting in a change of the quark content 
of particles in a medium. Consequently, the properties of the charge fluctuations seen 
in Figs. 4-6 are only necessary but not sufficient conditions to verify the existence of a 
critical endpoint (CEP) in the phase diagram. To verify the existence of the CEP one 
would need to observe a non-monotonic behavior of the net quark number (or electric 
charge) susceptibilities as shown in Fig. 7. The non-monotonic behavior of Xq along the 
boundary line appears since the quark number fluctuations are finite along the 2nd and 
the 1st order transition and diverge at the TCP. 

The finite structure of the charge susceptibilities along the first order line is valid 
under the assumption that this transition appears in equilibrium. Admitting deviations 
from the idealized equilibrium situation, the first order phase transition is intimately 
linked with the existence of a convex anomaly in thermodynamic pressure [20]. There 
is an interval of energy or baryon number density where the derivative of pressure 
dP/dV > 0. Such an anomalous behavior describes the region of instabilities bounded 
by the spinodal lines where the pressure derivatives taken at constant temperature or 
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Figure 6. The ratio of fourth- to second-order cumulants of quark number at /Zq = 
(left-hand figure) . The quark number density Uq with respect to the net quark number 
fluctuations Xq (right) as a function of T/Tq for various ^iqjT . Horizontal lines show 
the infinite temperature ideal gas values and the HRG model prediction (solid lines). 
The LGT results are from [H fT7] 



entropy vanishes. Changing the temperature or entropy results as the isothermal or 
isentropic spinodal line in the (T, ng)-plane with being the net quark number density. 
The spinodal decomposition of the system was argued to enhance the baryon number 
and strangeness fluctuations [20] . 

Fig. 8-left shows the phase diagram obtained in the NJL model formulated at finite 
current quark mass that accounts for spinodal instability. From this figure it is clear 
that the isothermal spinodal lines end up at the CEP. It is thus natural to explore 
how the charge fluctuations develop when going beyond the critical end point in the 
direction of the first order phase transition. Fig. 8-right shows the evolution of the net 
quark number fluctuations along the path of fixed T = 50 MeV in the (T, r;,g)-plane. 
The resulting behavior is quite interesting. When entering the mixed phase there is 
singularity in Xq that appears when crossing the isothermal spinodal lines, where the 
fluctuations diverge and change the sign. Between the spinodal lines the susceptibihty 
is negative. This implies a strong instability in the baryon number fluctuations when 
crossing the transition region between chirally symmetric and broken phase. 

The behavior of Xq seen in Fig. 8 is a direct consequence of the following 
thermodynamic relation 

'dP\ 

dVJr' VXq 

what connects the pressure derivative with the charge density Uq and the corresponding 
susceptibility Xq- Along the isothermal spinodal lines the pressure derivative in ([1]) 
vanishes. Thus, for non vanishing density rtq its fluctuations have to diverge to satisfy 
the relation ([T]). In addition, since the pressure derivative dP/dV\T changes sign when 
crossing the spinodal points, there should be discontinuity and the corresponding sign 
change in the divergence of Xq seen in Fig. 8. Due to linear relation between Xqi 



A =-^7-. (1) 



QCD phase diagram and charge fluctuations 



7 




Figure 7. The quark number (left) and isovector (right) susceptibiHties Xq Xi 
as functions of the temperature along the phase boundary. In the left-hand figure the 
solid (dashed) line denotes Xq the chirally broken (symmetric) phase. The vertical 
dotted-line indicates the position of the tricritical point TCP. The calculations were 
done in the chiral limit in the NJL model [9j. 



Xi and the electric charge susceptibihty xq a similar behavior as seen in Fig. 8 for Xq 
should be also there for xq- 

From Fig. 8 it is clear that when approaching the CEP from the side of the 
first order transition the region of instability shrinks and disappears completely at the 
CEP. Consequently, the CEP singularity in Xq appears from the matching of the two 
positive singular branches of Xq- Thus, the singular properties at the CEP are the 
remnant of divergent behavior of Xq at the first order chiral phase transition. From the 
above discussion it is clear that divergent charge fluctuations are not only attributed 
to the second order critical end point or TCP but are also there at the first order 
transition if the spinodal phase decomposition sets in. Consequently, the fluctuations of 
the conserved charges could be considered as a signal of the first order transition that 
is expected in the QCD phase diagram [5T]. 

In summary, we have presented a brief discussion of the phase structure expected 
in two flavor QCD. We have discussed the properties of charge density fluctuations 
when crossing the phase boundary. It was shown that for equilibrium transition the 
fluctuations of the net baryon number and electric charge exhibit a non-monotonic 
behavior along the QCD phase boundary that could be applied in heavy ion collisions 
as the signal for the existence of the critical end point. We have also argued that 
admitting deviations from the equilibrium picture for the first order transition the above 
non-monotonic structure is to be modified. In the presence of spinodal instabilities the 
fluctuations of the conserved charges diverge when crossing the isothermal spinodal lines 
of the 1st order transition. Consequently, such fluctuations can be used to probe the 
first order phase transition in heavy ion collisions. 
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Figure 8. The left-hand figure: The phase diagram in the temperature T and quark 
number density Uq plane in the NJL model. The filled point indicates the CEP. The full 
lines starting at the CEP represent boundary of the mixed phase in equilibrium. The 
broken-curves are the isothermal whereas the dashed-broken ones are the isentropic 
spinodal lines. The right-hand figure: the net quark number density fluctuations 
Xq/A^ normalized to momentum cutt-off A as a function of the quark number density 
Uq across the first order phase transition. The Xq was calculated in the NJL model 
along the line of constant temperature T = 50 MeV. 
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